General Relativistic Contributions in Transformation Optics 
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One potentially realistic specification for devices designed with transformation optics is that they 
operate with high precision in curved space-time, such as Earth orbit. This raises the question 
of what, if any, role does space-time curvature play in determining transformation media? Trans- 
formation optics has been based on a three-vector representation of Maxwell's equations in flat 
Minkowski space-time. I discuss a completely covariant, manifestly four-dimensional approach that 
enables transformations in arbitrary space-times, and demonstrate this approach for stable circular 
orbits in the spherically symmetric Schwarzschild geometry. Finally, I estimate the magnitude of 
curvature induced contributions to satellite-borne transformation media in Earth orbit and com- 
ment on the level of precision required for metamaterial fabrication before such contributions become 
important. 



^ • PACS numbers: 42.15.Eq, 42.70.-a, 78.67.Pt, 81.05.Xj, 91.10.Sp 



I. INTRODUCTION 



C/3 

^ ' Since its recent inception, metamaterial-based transformation optics has shown great promise as a new approach 
for manipulatingelectromagnetic fields and designing, novel optical and electromagnetic devices, such as the electro- 
magnetic cloak , diffractionless "superlens" 0, 0] , and temporal cloaking fs', y\ ; for a review, see El • The 
original approach to transformation optics used purely spatial transformations [1, 12.[l3j. A step toward generalizing 
transformation optics came from matching the constitutive relations for electromagnetic fields propagating in curved 
vacuum space-times with the constitutive relations for fields in stationary materials residing in Minkowski space-time 

[IMl. 

^ '' Common to these initial approaches is that they are based on a 3-vector representation of Maxwell's equations in 
O flat, Minkowski space-time. This excludes high precision applications that must take account of general relativistic 
curvature effects. Suppose, for example, that one would like to build an orbiting telescope using a superlens based 
_ 1 I on a space- folding transformation [14]. Would the well-known transformation media required for such a lens in 
^ ' flat space work equally well in orbit or does it require some modification to perform the same function? Another 
, example might be a very sensitive satellite antenna designed with transformation optics [l^ • The 3-vector approach to 
7-H ■ transformation optics cannot account for space-time effects, but a further gen eralization of transformation optics has 
lO ' recently been developed that extends the limits of previous approaches p^.ll9l|. This completely covariant, manifestly 
four-dimensional theory has already been shown to encompass all transformations in a unified framework that can 
^sj \, also accommodate relative motion IToj and arbitrary, non- vacuum, initial linear dielectric embeddings ^2B\. Here I 
. demonstrate that this approach also allows for transformations in curved space-times, thereby accounting for possible 
' general relativistic contributions to the resulting transformation media. 

. By using a curved space-time metric instead of the Minkowski metric it is possible to obtain an exact result, but 
^ ' because the scenario under consideration (Earth orbit) is decidedly non-relativistic it is advantageous to look at the 
non-relativistic limit. In particular, we could ask the question: If we design a device in fiat (or "Newtonian") space- 
time, how different does the transformation medium need to be for the device to actually work as desired in Earth 



' orbit? Given the weakness of gravity, we should expect the answer to be that the necessary material parameters are 
■ ■ precisely those expected from the fiat space-time calculations, plus something small, plus some higher order terms. 

To be more concrete, the envisaged scenario is that of a device in stable circular orbit around a spherically symmetric 
massive object whose space-time can be described by the Schwarzschild metric. We will restrict our attention to 
radially ingoing waves, and consider transformations over some small solid angle. This last restriction ensures that 
the device can be described in a single locally fiat patch by a co-moving observer. 

There is always a choice of coordinate system when specifying optical transformations. But in transformation 
optics, where one wants to specify a transformation that translates into a useful device, some coordinates are more 
appropriate than others, such as cylindrical coordinates for a cylindrical cloak Q. This will of course depend on the 
application in mind, but for the manipulation of radially ingoing light near a massive body, Schwarzschild coordinates 
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are the best choice. If one wanted to talk about transformations in the local frame of the device for fields known in the 
local frame of the device then everything is local and fiat, but does not describe the situation under consideration here. 
One could express the incoming fields in the local frame and then talk about everything locally as just mentioned. 
That approach is valid, but still requires all the same machinery developed here, and has the disadvantage of losing 
the intuitive description of optical transformations that make transformation optics so useful. 

There are some complications involved that make identifying and isolating the general relativistic effects more 
difficult. First is the question of how to calculate the purely Newtonian results with which to compare the full results. 
The purely Newtonian results would be those found for a device orbiting around a massive spherically symmetric 
body in flat space-time, i.e. in a theory of gravity where the space-time is not curved by the massive body. In flat 
space-time, the transformation medium for a device moving in rectilinear motion would be obtained through a Lorentz 
boost [l^. However, the orbiting device is non-inertial, being accelerated through its orbit, and the full orbital motion 
cannot be described by a boost. The best we can do is boost to an instantaneously co-moving frame at some point 
on the orbit. Thus, Sec. IV Al compares the expanded full results with the instantaneously co-moving device obtained 
via a boost in flat space-time. 

Secondly, the Newtonian limit of the material parameters will still depend on t he Ne wtonian velocity of the device, 
which for a circular orbit at radius r around a spherically symmetric mass M is Mjr. To clarify this point, consider 
for the moment a cloak moving with velocity v through a laboratory in flat space-time. The fields measured in the 
lab frame are not the same as the fields seen by the moving cloak. So one must either express the fields in the cloak's 
frame (via a boost) and do the transformation there, or do all the calculations in the lab frame and then boost the 
resulting material parameters to the cloak's frame. Expanding the final material parameters (in the cloak's frame) in 
the velocity v finds the same cloak parameters for a stationary cloak, plus some terms that are linear in plus some 
higher order terms [l9l |. In this example, the term linear in v could be interpreted as a velocity-induced correction to 
the stationary cloak. Returning to the scenario of an orbiting device, velocity dependent terms make up part of the 
Newtonian result, and care must be taken to isolate the curvature contributions. It should be emphasized that because 
the orbital velocity is expressed in terms of M, not every occurrence of M in the expanded material parameters is a 
curvature contribution. 

This article is organized as follows. Section |ll]briefiy summarizes the covariant, four-dimensional theoryof electro- 
dynamics. Section IlIII summarizes the covariant approach to transformation optics developed in Refs. (isl. [l9j. and 
presents the main result found therein as Eq. (|10|) . Section IIVI demonstrates how this approach to transformation 
optics may be used in a curved space-time; in particular, for a stable circular orbit of the Schwarzschild space-time. 
First I discuss how to express the permeability, permittivity, and magnetoelectric couplings in a local frame co-moving 
with the material, where these parameters take on their usual physical interpretation. This is then demonstrated with 
some representative transformations. Section |V] addresses the question of how these results are different from those in 
fiat space-time and shows that the results cannot be due solely to the orbital motion. In the non-relativistic, Newto- 
nian limit, the results agree with those of a boost in fiat space-time, as expected. Higher order terms of the expansion 
contribute curvature induced corrections. I conclude with Sec. |Vll I use the geometrized units where G = c = 1. 



II. CLASSICAL ELECTRODYNAMICS 

For details of covariant electromagnetic theory one may consult a number of excellent sources such as (2T| - [2^ , here 
I describe the salient features of the theory required for our purpose. Assume space-time to consist of a manifold M 
and metric g. The electric field E and magnetic fiux B are combined into a single mathematical object, the covariant 
field strength tensor F, that in a local Cartesian frame or Minkowski space-time has the matrix representation 

/ ~Ey —Ez\ 

p _ \ Ex Bz -By , 
-^M" ~ I _ R n R I • y'-) 

Additionally, the electric flux D and magnetic field H are combined into the covariant excitation tensor G, that in a 
local Cartesian frame or Minkowski space-time has the matrix representation 
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Maxwell's equations are succinctly expressed as dF = 0, and dG = J, where d is the exterior derivative, and J is the 
charge-current 3- form [2l| . 
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In a linear dielectric medium there exists a relationship between F and G given by a constitutive relation that may 
be expressed as [H, [l^ 

G = x(*F), (3) 

which in component form reads 

Gfiy = X^zy"'^ *a/3'^'' ^^rp- (4) 

In Eq. ([2]), * is the Hodge dual on (Af, g), which for present purposes is to be understood as a map from 2-forms to 
2-forms having components 

*a/'' = ^^e„;3,p5'^v^ (5) 

where ea/Sa-p is the totally antisymmetric symbol, and g is the determinant of the metric tensor Qf^^. Also note that 
g'^'^ is the matrix inverse of g^^jy. The tensor x contains information on the dielectric material's properties such as 
permittivity, permeability, and magnetoelectric couplings. We take x to be independently antisymmetric on its first 
two and last two indices, and in vacuum x(*F) = *F. This last condition means that the classical vacuum is treated 
as a linear dielectric with trivial x, recovering the usual, trivial, constitutive relations in vacuum. 

The components of the constitutive equation provide a set of six independent equations that in a local frame or 
Minkowski space-time can be collected in the form 

Ha = irX'B, + UnaEb, Da - inaEb + {l2*)a'B,, (6) 

where the notation a denotes a 3 x 3 matrix. Rearranging these to 

gives a representation that may be more familiar and in which subsequent results will be expressed. These three- 
dimensional representations of the completely covariant Eq. ^ are essentially equivalent, and are related by 

£ = £ - 72 Ai7i , 7i = -/^7i , 72 = 72 (8) 

However, one should be aware that these 3x3 matrices are not strictly tensors but sim ply components of x that have 
been collected into matrices. The components of x have been identified elsewhere Il9l|. 



III. TRANSFORMATION OPTICS 



To imderstand transformation optics, start with an initial space-time manifold (M, g,*), field configuration 
(F, G, J), and material distribution x, where dF = 0, dG — J, and G = x(*F). Imagine now a map T : M ^ M Q M 
that maps M to some image M and transforms the electromagnetic fields to a new configuration (F, G, J). Because 
the underlying space-time is physically unaltered the manifold is still described by (A/, g,*). But for the new field 
configuration to be physically supported, there must exist a new material distribution x- Therefore dF = 0, dG — J, 
and G = x(*F) holds on M. Such a transformation could, for example, map M to an image M that contains a hole, 
i.e. a region from which the fields will be excluded, as in the case of an electromagnetic cloak. 

Using the inverse, T, of the map T, the initial F and G are related to the final F and G by the puUback of T, 
denoted as T*. This implies 



G = r*(G) = r*(x(*F)) = x(*r*(F)), 

which can be solved for % as a function of a; G M, giving [l^, [l^ 



up 



T(x) 



r(x) 



(9) 



(10) 



In Eq. ([TU| A is the Jacobian matrix of T, A~^ is the matrix inverse of A, and in solving for x we have made use of 
the fact that on a four-dimensional Lorentzian manifold, acting twice with ★ returns the negative, F = — F. The 
initial material tensor x ^-nd first * must be evaluated at Tix), while everything else is evaluated at x. Equation (jlOp 
represents the core of transformation optics in linear dielectric materials. 

Note that Xa^" on the right hand side of Eq. pU)) need not describe the vacuum, so this approach can be used 
for transformation optics in arbitrary, non- vacuum, initial linear dielectric media [20| . Furthermore, the space-time is 
unspecified. Given a desired space-time, such as near Earth, the appropriate metric enters Eq. (|10l) through 7k-, given 
by Eq. ([S|), as is demonstrated below. 
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IV. TRANSFORMATION OPTICS IN CURVED SPACE-TIME 



To demonstrate how Eq. (|10p may be applied in a curved space-time, consider a situation where the desired device 
travels in a stable circular orbit around a spherically symmetric massive body, such as a planet. Although in this 
example the planet is not rotating, it illustrates the idea of a satellite-borne device operating in Earth orbit. The space- 
time around the planet is described by the Schwarzschild metric, which in the Schwarzschild coordinates (t, r, 9, if) 
takes the matrix representation 



/-(l-^) \ 

(1-^)"' 

r2 

V r^sin^e/ 



(11) 



where M is the planetary mass. The contravariant version of the metric, 5'^'^, as needed in Eq. ([5]), is the matrix 
inverse of 17^1,. 

The utility of transformation optics rests on the specification of a transformation that translates into a device, 
and one is always free to choose the coordinate system in which to specify a transformation. However, certain 
coordinates lend themselves more readily to particular situations, such as cylindrical coordinates for designing a 
cylindrical cloak. The physical system I have in mind is of light coming in from space. The simplest thing to look at 
is radially ingoing light with frequency defined with respect to Schwarzschild time t. I want the radially ingoing light 
to be directed or behave in a particular way, from the Schwarzschild coordinate perspective, that can be specified 
through a transformation. Thus it is assumed that the most appropriate coordinates to work in are the Schwarzschild 
coordinates. 

Geometrizing the material parameters into the tensor x means that x can be used in a general tensor framework, 
but interpreting the components of x permeability, permittivity, etc. is understood, measured, and built, in a flat, 
local frame (or Minkowski space-time) according to Eq. ([7]). So while we might want to work in the Schwarzschild 
coordinates and define the optical transformation there, these coordinates don't make any sense for talking about 
dielectric materials if you actually wanted to build something, thus we have to express the transformation media in 
the locally flat frame of an observer orbiting with the device. The orbiting observer has his own notion of time with 
respect to which he is defining spatial hypersurfaces; essentially determining his own breakdown of the field strength 
tensor and his notion of material. In such a locally flat frame the metric becomes Minkowskian, and the components 
of X may be given an interpretation in terms of the familiar material properties of Eqs. An observer co- moving 
with this frame would measure and interpret the material properties in the same way as in a laboratory on Earth. 

A. Construction of a local frame 

In flat space-time, one may readily switch between inertial observers by using a Lorentz boost, but switching 
between inertial observers in a curved space-time is more complicated. To express x in a local frame co-moving with 
the device, we must first specify the motion of the local frame in Schwarzschild coordinates and then construct a 
boost-like operation that allows us to switch between the global Schwarzschild coordinates and the local frame. Thus, 
the local frame is essentially the set of Cartesian axes carried by the co-moving observer. 

For simplicity, and without great loss of generality, we may restrict attention to motion in the equatorial plane, 
= 7^/2. The four- velocity of a small object in a circular orbit of radius r, in the Schwarzschild coordinates, is [1^ 

(12) 

'1 "I 

Note that in the Newtonian limit r ^ M this reduces to the usual Newtonian angular velocity 

I [m 

iON = -\ — (13) 
r \ r 

of a particle in circular motion about a mass M. Stable timelike circular orbits only exist for r > 6Af, so the 
potentially problematic value r — 3M is not under consideration. 

Let Co be a proper-time-directed basis vector in the local frame ca- The capital Latin indices of the local frame 
are raised and lowered by the Minkowski metric ?7ab- In the local frame, the worldline of the co- moving observer 
is directed along the direction of proper time. It follows that in Schwarzschild coordinates, Bq = u'^. We must now 
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find spatial unit vectors ca such that e^e^ = tjab- Because only Cq is constrained, there is some choice as to the 
orientation of the spatial frame, and the following choice of frame will be discussed further in Sec. IV CI Let the local 
observer's x-axis coincide with the r-direction at all times, which implies 




A suitable choice for 62 is 



0,0,i,0 



(14) 



(15) 



Because ei has only an r-component and 62 has only a ^-component, 63 can only have t and ip components. This, 
together with the orthogonality conditions, determines 



1 



Mr y/r ~ 2M 



(16) 



^ v/(r - 3 A/) \^ V ~ 2^^) ' / 

A transformation matrix for switching between the global Schwarzschild coordinates and the local frame may be 
instructed with the basis vectors so determim 
UA in the local frame by ua — S^'^n^, where 



constructed with the basis vectors so determined. A 1-form in Schwarzschild coordinates is transformed to a 1-form 



I _ 3M 







M 



-2M 



1_ 3M 
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-Lji-m 



(17) 



while vectors are transformed by S^^, the transpose of the inverse of S'^''. 

One can readily verify that these matrices transform 1-forms and vectors as desired. In particular, transforming 
the Schwarzschild metric to the local frame by 



VAB = 



0\ 
10 
10 
0001/ 



(18) 



returns the Cartesian Minkowski metric, as required (in the equatorial plane). I emphasize that in this article we 
are dealing with two different operations that, perhaps due to an inadequacy of language, are both referred to as 
"transformations" . In this subsection we have constructed the matrix operation Eq. (jl7|) that allows us to transform 
tensors between the global Schwarzschild coordinates and the local co-moving frame of the orbiting device, i.e. taking 
a physical object, such as a vector, expressed in Schwarzschild coordinates and re-expressing it in the local coordinates 
of the co-moving observer. This frame transformation is distinct from the map T of transformation optics (see Eq. 
(|ni)) that is typically given in terms of a coordinate transformation. I will subsequently distinguish these as "frame 
transformation" and "optical transformation" . 



B. Optical Transformations in Schwarzschild Space-Time 



A local, flat, Cartesian frame has been constructed in which the transformation medium parameters have meaning 
for a CO- moving observer, in the sense of Eq. ([T]). As is evident from Eq. ([9]), the optical transformation that must be 
specified is a map T : M C M — >■ M, which is a transformation on the Schwarzschild coordinates of the form 

T{t, r, cp) = it'., /, 0', cp') = (/o(i, r, 9, ip), /i(t, r, 0, cp),f2it, r, 6, ip), hit, r, 9, ip)) . (19) 

Calculating the Jacobian matrix of T and using Eq. pop leads to an expression for x, which may be frame transformed 
to the local frame by 

Xab'"^ — ^ A^ ^ b" pX^iJ^^ ■ (20) 
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The material parameters may be extracted from x and expressed in the representation of Eq. ([T])- Note that while 
X on the right hand side of Eq. (^0]) is expressed in Schwarzschild coordinates {t,r,9,ip), x on the left hand side 
is expressed in the local Cartesian coordinates (x'^ , , x'^ , x^) = {T,x,y,z). Generic optical transformations such as 



Eq. (jl9p are computationally intensive, and yield complicated results. Consider instead some simpler representative 
optical transformations. 



1. Radial Optical Transformations 



One application of near Earth transformation optics that springs to mind (disregarding any issues related to near- 
field or far-field applicability, or boundary conditions) is that of an orbiting telescope fashioned from a diffractionless 
perfect lens, or "superlens" @, 0|- A perfect lens of width d, in the Cartesian coordinates of flat space-time, has been 
frequently discussed in terms of a linear space-folding optical transformation [26l - l28j such as 



T{t,x,y,z) = {t',x',y',z') = 



{t,x,y,z) 
{t,-ax,y,z) 
{t, X — 2ad, y, z) 



a; < 0, 

< a; < d, a > 0, 
X > d, 



(21) 



Extrapolating to Schwarzschild space-time, one could imagine a perfect lens in Earth orbit designed to focus radially 
ingoing light rays, over some small solid angle, with a radius-folding optical transformation. However, it has been 
argued that because the space-folding transformation does not preserve orientation, certain subtleties in the TO 
procedure mean that this transformation does not adequately describe a superlens [29] . Consider instead the general 
class of orientation preserving radial optical transformations of the form 



T{t,r,e,^)^{t',r',e\^')^{t,f{r),e,^). 



(22) 



such that the resulting transformation medium occupies some appropriate radial region and small solid angle. The 
Jacobian matrix of the transformation is 



A ^ - 



/I ON 

/,r 

10 

VO ly 



(23) 



and ^^^^'^ is obtained for the Schwarzschild space-time by plugging the Schwarzschild metric, Eq. (jlip . into Eq. ([5]). 
The Jacobian matrix A"^ and ^^^^ are now used in Eq. (|10|) to calculate the transformation medium X\^^ ■ Note 
that because the optical transformation takes place in vacuum, the initial XaB^" i the right hand side of Eq. (fTO|) , 



is simply the vacuum x^ so that Eq. pOl) reduces to 



crp 



T[x) 



(A-i)^(A-i)^ 



(24) 



Next, this calculated value of X\J^ is used together with Sj^^ in Eq. (|20|) to find Xab~^^ ^ which is the material 
parameter tensor x expressed in the local, co-moving Cartesian frame. 

Once Xab^^ known, it would be nice to give the material parameters their usual interpretations in terms of 
permeability, permittivity, etc. The correspondence between the components of x^p^^ and the components of the 
matrices e*, fi^^ ^ 7]*, and 72 in the representation of Eq. ([6]) are given explicitly in 18], or in the Appendix of [l^. 
Once £*, 7i, and 73 have been identified, they are then converted, through application of Eqs. ([5]), to the usual 
e, fi, 71 and 72 of the representation given by Eq. ([T]). These are finally the material parameters that realize the 
optical transformation Eq. (|22|) . Expressed in the local, co-moving Cartesian frame, they are 




fr{r - 3M) 
r3/ - M(2r3 -|- p) 
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(25) 



7i 




(26) 
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with fi = £ and 72 =72- This shows that space-time curvature effects appear even for the relatively simple radial 
transformation of Eq. (22). In this case it is straightforward to see that in the limit r ^ M and for localized 
transformations /(r) ^ r this reduces to the Newtonian result 



'l/f 0~ 
£ - I /' 
/' 



71 =0. 



(27) 



Typically, dielectric motion requires magnetoelectric coupling terms and also tends to change the required values 
of £ and ^ [ll, 1131 , indeed this motional effect has been used to interpret magnetoelectric couplings as velocities in 
transformation optics [14]. But this does not happen for the optical transformation Eq. ((22|) . even though the device 
is in relative motion throughout its orbit. Why not? Consider for a moment electromagnetic fields measured in a lab- 
oratory in flat space-time. If a dielectric medium moves through the lab with speed then the fields experienced by 
the medium are not those measured in the lab frame, but are instead the lab-measured fields boosted to the medium 
frame. By boosting to the medium frame, one finds that medium traveling perpendicularly to the transformation 
direction requires no modification over stationary medium, while medium moving parallel to the transformation direc- 
tion requires velocity dependent modifications to the material parameters of the stationary medium. Returning now 
to the orbiting device, we have chosen an orbiting frame such that the instantaneous velocity is always perpendicular 
to the r-direction, the direction of the action of the optical transformation. Therefore, in this particular case we 
should not expect any contribution from the velocity. 

Mass dependent, general relativistic contributions start to appear even for the relatively simple transformation 
Eq. (22). As a slightly more general class of radial transformations, consider 



(28) 



such that the equivalent material occupies some appropriate radial region and small solid angle. In this case the 
corresponding material parameters, in the local Cartesian frame and occupying the appropriate spatial region, are 

. _ <7(g-2M)(r-3Af) 
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(2M-r)g,, 
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(2M-r)[Af(2r'+g')-gr^]g,^ 
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(29) 
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-T 
72 



v/M(^^^2M) 



where everything is evaluated on the equator, 

H — gr{r^ — 

and 



H 








3M)g^g,rg,v 




2M 



2M 



(2M(r3 - g^)) 



2 2 
9 5,v>, 



(30) 



(31) 



A^{g- 2M) [r'^g - M{2r^ 



Mg'gl 



(32) 



(33) 



where everything is evaluated on the equator, and 

A^r[r- 2M)[r - 3M) - Mg^^. 

Due to the presence of M in Eqs. (|29l) and (150]) it is clear that there is some additional contribution over what might 
be expected from a similar optical transformation leading to a stationary device in flat space-time. However, distin- 
guishing the contribution due to the orbital velocity from that of the space-time curvature requires some additional 
effort. Before turning to such an analysis, consider an optical transformation that mixes space and time components. 



2. Optical Transformations Mixing Space and Time 

Recently, optical transformations mixing space and time have generated some interest [l9l Isij , and have been 
examined in the context of frequency converting active metamaterials [3l| . The transformation discussed in Ref. [sij 
was given by 



Tit, X, y, z) = {t', x\ y', z') = ( — ^—r,x, y, z] 

\ax -\- J 



(34) 
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within a stationary slab in flat space-time. For our purposes, consider a similar optical transformation, in the 
Schwarzschild coordinates. 

Tit, r, 9, ip) = {t\ r', 6', = (/o(t, r), r, 0, ^) , (35) 

piecewise defined over a suitable radial region and small solid angle. It is assumed that the portions of the optical 
transformation outside this small active region is boundary matched and leaves the external vacuum region unaffected, 
but may not be smooth at the boundary. In this case the corresponding material parameters, in the local, co-moving, 
Cartesian frame, are found to be 



(3M - r)/o,t 








(1- 



2M^ 



(3M - r)/o,t 
v/M(r - 3A/)/o,. 



(1 



2M\ 







3M)/o 



B ~ M (1 



2M'\ 




(36) 
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72 



Vr - 2M 
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where 
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(1 - 2M) V^—3Mfo,thr 



B = M-{r- 2M)fl 







1 - 2M) ^/V^iMh.th, 




(37) 



(38) 



This demonstrates another example where a class of optical transformations in an orbit of Schwarzschild space-time 
requires an M-dependent transformation medium. Again the question remains of how to disentangle the velocity and 
curvature contributions, since both of these will appear as M dependencies. More precisely, how can one determine 
the curvature induced contribution over what one would obtain for a device performing the same action when in orbit 
around a Newtonian planet of mass M residing in a fiat, Minkowskian background? We will return to this question 
in Sec. |Vl but first consider one more example of interest. 



3. Orbiting Cloaking Device 

Suppose it is desired to construct an orbiting cloaking device that hides a region from radially ingoing light rays. 
Let the device occupy a segment of an annulus that circumscribes a sphere of radius Tq, and let the interior region 
be of a size to contain a sphere of radius rj, as in Fig. [T] Let coordinates on the vacuum manifold be {t' ,r' ,9' ,(p') 
and coordinates on the manifold containing the cloak be (t, r, 9, ip). For simplicity we may consider a two-dimensional 
cross-section of the cloak lying in the equatorial plane 9 = 7r/2. The three-dimensional cloak can obtained by either 
by rotating the two-dimensional cross-section about the radial line passing through its center, or, for a different cloak 
geometry, rotating the two-dimensional cross section through an azimuthal angle A9, such that the inscribed circles 
of Fig. [T] become arced cylinders. Let the center of the cloak orbit at a radius Tc, and be instantaneously located at 
</j = 0. A radially ingoing ray enters the cloak at a point ri = Vc + Tq and exits the cloak at point r2 = Vc — Tq, as in 

Fig.m 

Instead of following the chord connecting ri and r2, the ray is deflected upon entering the cloak and follows a 
diverted path. Such a path may be described via a map T : M ^ M from the chord in the vacuum manifold to the 
manifold containing the cloak. This may be thought of as a coordinate transformation that rotates a point (r',ip'), 
lying on the chord, to a point (r ~ r' ,ip{r' ,(p')) on the diverted path. The maximum amount of deflection from a 
given chord will depend on the ingoing angle Lp' , as in Fig. [51 An ingoing ray with ip' — Q has a maximum deflection 
equal to the half-angle spanned by the inner spherical region (e.g. small dotted circle in Fig. [T]) 

(39) 

An ingoing ray with p' — uJo, where 




FIG. 1. A radially ingoing ray enters the cloak at the point r\ = + To and exits at the point r2 — Vc — ro- Instead of following 
the chord connecting ri and r2 (dashed line), light is deflected upon entering the cloak and follows a diverted path (solid line). 
The diverted path is described by a bump function and is therefore C°° at the cloak boundary. 




FIG. 2. Radially ingoing light rays follow a deflected trajectory as they pass through the cloak. The cloaked region (shaded) 
is not spherical, but is designed to contain a sphere of radius ri and angular size 2(jJi. The deflection amplitude varies with the 
angle of the ingoing ray and has maximum value tOi for = and decreases linearly to zero at ip — ujo. 



is the half-angular size of the cloak (or the sphere circumscribed by the cloak - e.g. large dotted circle in Fig. [TJ, 
grazes the exterior surface of the cloak and will not be deflected at all. 
To find a map satisfying these conditions, let (p be given by 



ip{r' = r, ip') = tp' + iuj. 



where 



b{r) = 



cxp 




1 - 



(»'c - To) <r < (re + To), 
elsewhere 



(41) 



(42) 



is a "bump function" of unit amplitude. Although a bump function is piecewise defined, it is C°° everywhere and 
therefore smooth across the boundary. For calculational purposes assume ip' > 0, allowing the cloak parameters to 
be specified in the upper part of Fig. [2] The lower part may be obtained by reflection. 

Recall from Sec. Illll that. while it may be convenient to describe the desired trajectory with a map T : A/ — > M, 
such as Eq. (j4ip . from the vacuum manifold to the manifold with material, it is the pullback of the inverse map 
T : A/ — > M that maps the electromagnetic field strength and excitation tensors F and G from the vacuum manifold 
to the manifold with material. It is only then that the corresponding trajectory of radially ingoing light rays passing 
through the cloak material will be diverted as in Fig. [51 Therefore, let 



r{t, r, e, ip) = {f, r', e', ^') = {t, r, e, Mr, ip)) 



(43) 
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where 

/3(^^) = f^S^- ^>-.6W. (44) 

From Eqs. (flO)) . (|17p and pO|) . the material parameters for the cfoak, expressed in the local Cartesian frame and 
representation of Eq. ([7]), are 







(2M - r) Vr(r - 3M)/3,. 





(3M - r)/3, 








C-r(2M-r)V| 



(45) 



where 



72 = 7i 



C 



-(1-/3%) 





-v/r(r-3M)/3,^/3,. 
^r(r - 3M)/3,^/3,. 



C = 2M-r + M/32^. 



(46) 



(47) 



Consider a specific example of such a cloak. The orbital radius of a typical satellite is Tc ~ lO^cm, and the 
geometrized mass of Earth is Me ~ 0.3cto. Let — lOOcm and Tq — SOOcto. The corresponding material parameters 
for such a cloak are plotted in Fig. [3l 






-O.OOOOll 
-0.000021 
-0.00003; 



(d) l^y 



(e) 7zH 



FIG. 3. Material parameters for half of a two-dimensional slice of an orbiting cloaking device described by the transformation 
given in Eq. H43|) . A three-dimensional cloak can be constructed in two ways: either by rotating about the r axis, or by reflection 
across the r axis and then stacking copies of the slice to form a tube-like hidden region. 



Because the angular size of the cloak is so small, the curved shape, greatly exaggerated in Figs. [T] and [H is not 
perceptible. The density plots of Fig. [3] look very much as they would for a similarly defined cloak in flat space, 
which returns us to the question posed at the end of the Sec. IIVB 2[ namely, how can one determine the curvature 
induced contribution to the material over what one would obtain for a device performing the same action in a flat, 
Minkowskian background? This is addressed in the next section. 
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V. DISCUSSION 



It was mentioned earlier that for transformation optics, a dielectric moving with respect to the frame in which 
the optical transformation has been defined typically requires magnetoelectric couplings, and that in cases where an 
optical transformation results in such couplings, they can sometimes be interpreted as a velocity [3, [l^. Because 
the calculations of Sec. IIVI have been demonstrated for a material in orbital motion, and because the four- velocity is 
expressed in terms of the planetary mass and orbital radius, it is justified to ask whether the Af-dependent results 
truly include space-time curvature effects or are merely a consequence of the dielectric's orbital motion, after all, the 
Riemann curvature tensor does not appear explicitly in any of the calculations. But the Riemann tensor need not 
appear explicitly to find effects that only appear in curved space-times. 

The ideal approach would be to compare with results obtained in the absence of general relativity, i.e. the transfor- 
mation media required for the same optical transformations applied in conjunction with a Newtonian theory of gravity. 
For this purpose, consider a "Newtonian" system consisting of the dielectric device in circular orbit around a spherical 
mass ill, where the background space-time is Minkowskian, i.e. a system with vanishing space-time curvature. This 
assumption combines the relativistic covariance of electrodynamics with Newtonian gravity, and essentially consists 
of a dielectric device moving in a circle of radius r in Minkowski space-time, with angular velocity given by Eq. (jl3p . 
We may still consider optical transformations in spherical coordinates, and must still express the dielectric material 
properties in the local frame of the device. 

At first glance it might appear that since the correct Newtonian velocity is obtained by taking the Newtonian limit 
r ^ M of Eq. (fT2|) . one should proceed by using the Newtonian limit of Eq. ([T7|) to perform the same calculations. 
This approach is incorrect. In flat space-time, frame transformations between inertial observers is accomplished with 
a Lorentz boost, but the Newtonian limit of of Eq. pT)) is not a Lorentz boost. A circular orbit in flat space-time 
is non-inertial, and for such a motion there is no clearly accepted method for frame transforming between observers 
over the entire orbital path. The best we can do is use a Lorentz boost to frame transform to an instantaneously 
co-moving frame. 



A. Comparison with Boosted Frames 



The first thing we can do is compare the fully relativistic results with those expected in a boosted frame in flat 
space-time, which will reveal whether or not the full results arise solely from the orbital motion. Since it is assumed 
that the transformations in the Schwarzschild coordinates occur over a very small solid angle, and the final results are 
expressed in a locally flat Cartesian frame, we may compare with an equivalent Cartesian transformation in Minkowski 
space-time. For definiteness, return to the transformation of Eq. ([35|) . Due to the orientation of the local Cartesian 
frame chosen in Eqs. (|14p - this corresponds to a transformation of the form 



r(i, X, y, z) = {t', x', y', z') = (/o(i, x), x, y, z) , (48) 
boosted in the z direction with some speed /3. A Lorentz boost in the z direction has matrix representation (2]| 



/ 7 7/3N 

10 

10 

\lP 7,, 



(49) 



where 7 = (1 — /3^)~^/^. This frame-transforms components of a 1-form (or lower tensor indices), Wa = w^i'^a, such 
that the frame with Latin indices moves with velocity /3z in the frame with Greek indices. Components of a vector 
(or upper tensor indices) are frame-transformed by — L°',,v^, where L"^ is obtained by setting /3 — > — /? in L^^. 

Equation pO|) is perfectly compatible with such boosts _1!}|. Applying Eq. pO]) and a boost to the transformation 
of Eq. (j48p in Minkowski space-time results in the material parameters 



(50) 



7i = 72 



/o,i 



~/3(l - k,t) 

2 ,/3(i-/o,0 -V^~Wh^r 



(51) 



^T^/O,: 
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By comparing Eqs. (|50l) and ([5T|) with Eqs. (|36|) and (|37)) . it is clear that the fully relativistic results Eqs. (p6)) and 
([57)1 cannot arise from a pure boost in flat space-time. Choosing 

P = ±y/M/{r -2M) (52) 

will recover the desired values of fixx, fj-yy, and jixy (for the negative root of /3) but not the other non-zero components. 
Choosing instead 

JAlVPiM -r)+ r^{r - 2M)fl^ 

will recover, if real, the desired value of jlzz but none of the other non-zero components. However, there is no boost 
speed that will recover the desired fixz or ^lyz components. This demonstrates that the fully relativistic results 
obtained in Eqs. p6| and (|37l) are not equivalent to a boost in flat Minkowski space-time, and are therefore not solely 
a consequence of the orbital motion - there are space-time curvature contributions. 



B. Newtonian Limit 



Section [V Al demonstrates that the fully relativistic results are, in general, not simply a result of the orbital motion 
but also include curvature effects. But the curvature contribution was not isolated. If scientists of the future wanted to 
send a transformation optics designed probe to a black hole, general relativistic effects would surely become important. 
What about for near Earth applications? In the geometrized units used here, mass may be measured in distance. 
The geometrized mass of Earth is roughly 0.3cm, while the radius of a typical satellite orbit is on the order of lO^cm. 
Additionally, the velocity of a satellite in orbit is much less than the speed of light, so /3 << 1. Thus the next thing 
we can do is look at an expansion of the results in the Newtonian limit. 

From Eq. (|52l) it should be expected that in the non-relativistic, Newtonian limits /3 — )■ and r » M (equivalently 
M — >■ 0), Eqs. (PSI) and (1571) are equivalent to the flat space-time results boosted with the Newtonian speed for circular 
orbits, /3 — ^ Mjr. Expanding Eqs. (|36|) and ()37|) to lowest order in M, and Eqs. (f50| and (fSTj) to second order in /3 
shows this is indeed the case. Continuing with the expansion, one finds curvature terms appearing in the next order 
of the expansion, in terms proportional to 

M [M 

— \/ — • 54 

r V r 

Therefore, it is found that the requisite material parameters, described in a local, co-moving, Cartesian frame can be 
described as 

e = £0 ei -j- £2 + O ( — j (55) 

where Eq are the material parameters for an equivalent transformation in flat Minkowski space-time, ei are terms 
proportional M/r, and £2 are terms proportional to Eq. (l54l) . While only the velocity contributes to £1, both velocity 
and curvature contribute terms in £2. 

Near the horizon of a black hole, these corrections can approach £1 « 0.5, £2 ~ .35, which can represent a significant 
contribution for transformation media with £ ^ 1. On the other hand, for a satellite orbiting Earth with r ~ lO^cm 
these corrections are on the order of £1 ~ 10^^ and £2 ^ 10^ With light, differences of 1 part in 10^"^ are not 
difficult to measure with an interferometer, but such tiny corrections are not likely to have great significance for most 
applications of transformation optics in Earth orbit. On the other hand, in gradient indexed materials comprising 
holes or voids in glass the void volume scales with the index of refraction, thus one would have to control the void 
volume to these levels. Such precision is certainly challenging but is not inconceivable, thus it may be that such 
corrections will be required for highly sensitive applications in the future. 



C. Frame Choice 



The particular frame described in Sec. IIV Al is not the only choice of local frame. In fact, since the local x-axis is 
always radially aligned, the chosen frame is actually spinning with respect to the distant stars. It may be of interest 
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to instead consider a non-spinning local frame. This kind of frame is said to be "Fermi- Walker" transported [2i|, [33] ■ 
Since in the Fermi- Walker case the x-axis will not always be pointed radially outward, the parametrization of the 
local frame will be time dependent and the required frame transformations will also be much more complicated. The 
spinning frame of Sec. IIV was chosen to simplify the demonstration of general relativistic corrections. 

It was shown in Sec. IV Al that the fully relativistic results are not, in general, a consequence of the translational 
dielectric motion alone, but could they simply be a gyroscopic effect of the spinning frame? Indeed, the Lorentz 
transformations consist of both boosts and rotations, and the fact that pure boosts do not form a closed subgroup of 
the Lorentz group gives rise to, for example, Thomas Precession [33,, jB4j. However, it is easy to see that Eqs. and 
(|37)) cannot be obtained as a combined boost and rotation in flat space-time. Such a result would be obtained in a 
similar manner to Eqs. (15(11) and ([5T|) . although in this case the Lorentz transformation is slightly more complicated. 

A general Lorentz transformation is of the form 

L = exp[a; ■ S + i ■ K] (56) 

where the 4x4 matrices Si and Ki are the generators of rotations and boosts, u; is a 3- vector denoting an angle 
about an axis of rotation, and 

^ = tanh"\/3)/3 (57) 



is the boost 3- vector 3j| . In this case we may restrict ourselves to a boost in the z direction and a rotation in the xz 
plane. Taking the non-relativistic, Newtonian limits shows that the limiting behaviour of Eqs. ([55]) and ([57]) cannot 
be recovered for any non-zero uj. Thus the full result is not simply an artifact of boosts and rotations, but does in 
fact include curvature contributions. 



D. Geometric Optics 



While I have calculated the general relativistic, curvature induced contributions to the material parameters, one 
might also like to compare the difference in the propagation of light through the corrected and uncorrected materials 
by solving Maxwell's equations through the medium, e.g. by ray tracing through the cloak of Sec. IIV B 31 Of course, 
for the full material parameters of Eqs. (|45l) and (|46|) . the rays follow the trajectories of Fig. [T] hy design, but the 
question for a TO designed instrument would be whether the uncorrected parameters are "good enough" for the desired 
application. This means it would be desirable to compare the exact result with the ray trajectories through a device 
constructed from the material parameters obtained in the Newtonian limit, but still orbiting in the Schwarzschild 
geometry. Since GR contributions are generally weak, the material parameters and boundaries generally complicated, 
and the space-time curved, the calculation would likely require a more complicated dislocation approach f35'| rather 
than an eikonal or plane-wave ansatz in order to retain sufficient sensitivity in the geometric optics limit. Such a 
calculation lies outside the scope of the current article, the purpose of which is as the first foray of transformation 
optics in curved space-times; I therefore leave this for further study. 



VI. CONCLUSIONS 

The completely covariant approach to transformation optics introduced in Refs. [3, [l^ allows for arbitrary back- 
ground space-times. This enables space-time curvature effects to be accounted for in transformation optics. I have 
demonstrated the particular case of optical transformations for transformation media in a stable circular orbit of the 
spherically symmetric Schwarzschild geometry, simulating an Earth-orbiting satellite. In order to have their usual 
physical interpretations, the material parameters must be expressed in a locally flat Cartesian frame, co-moving with 
the satellite. After constructing such a frame, and the matrices used to frame transform % into it, some generalized 
optical transformations were considered to illustrate the calculations. As a specific example, a unidirectional orbiting 
cloak that redirects radially ingoing rays around a central cavity was considered in Sec. IIVB 31 

Of particular concern is whether the mass contributions to the results are a consequence of the space-time curvature 
or the orbital velocity, as both depend on M . However, it was explicitly demonstrated for the optical transformation 
given by Eq. psp that the result is not equivalent to motion without curvature, i.e. motion in flat space-time. The 
curvature effects were further disentangled by considering the non-relativistic, Newtonian limits. In these limits, the 
results are in fact obtainable as a pure boost in flat space-time, with no rotation. Higher order terms of the expansion 
provide curvature induced corrections. 

The transformation considered in detail, Eq. psp . generalizes a transformation that has been previously studied in 
the context of frequency shifting metamaterials. However, a similar frequency shifting effect, called the gravitational 
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redshift, arises naturally from the gravitational potential [2lj and is easily measurable on Earth [36[. It is, therefore, of 
no surprise that when overlaying a frequency altering transformation on the naturally occurring gravitational redshift, 
space-time curvature effects will be required to end up with the correct, desired frequency. Such corrections are small, 
but are potentially meaningful and accessible for future high precision applications of transformation optics. 
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